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Abstract. We apply results of |BP07I IBBBPOS] to compute limiting prob- 
ability profiles for various quantum walks in one and two dimensions. Using 
analytical machinery we show some features of the limit distribution that are 
not evident in an empirical intensity plot of the time 10,000 distribution. Some 
conjecutres are stated and computational techniques are discussed as well. 



1. Introduction 

The quantum walk on the integer lattice is a quantum analogue of the discrete- 
time finite-range random walk (hence the abbreviation QRW) . The process was first 
constructed in the 1990's by |ADZ93| . with the idea of using such a process for 
quantum computing. A mathematical analysis of one particular one-dimensional 
QRW, called the Hadamard QRW, was put forward in 2001 by ABN+OT. Those 
interested in a survey of the present state of knowledge may wish to consult [KemOS] 
as well as the more recent expository works |K en07l IVA081 IKonOSj . Among other 
properties, they showed that the motion of the quantum walker is ballistic: at time 
n, the location of the particle is typically found at distance 9{n) from the origin. 
This contrasts with the diffusive behavior of the classical random walk, which is 
found at distance 9{^/n) from the origin. A rigorous and more comprehensive anal- 
ysis via several methodologies was given by |CIR03| , and a thorough study of the 
general one-dimensional QRW with two chiralities appears in [BP07j . A number 
of papers on the subject of the quantum walk appear in the physics literature in 
the early 2000 's. 

Studies of lattice quantum walks in more than one dimension are less numerous. 
The first mathematical such study, of which we are aware, is IIKK04] , though 
some numerical results are found in MBSS02 . Ballistic behavior is established 
in [IKK04) , along with the possibility of bound states. Further aspects of the 
limiting distribution are discussed in WKKK08]. A rigorous treatment of the 
general lattice QRW may be found in the preprint BBBP08J. In particular, 
asymptotic formulae are given for the n-step transition amplitudes. Drawing on 
this work, the present paper examines a number of examples of QRWs in one and 
two dimensions. We prove the existence of phenomena new to the QRW literature 
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as well as resolving some computational issues arising in the application of results 
from [BBBP08] to specific quantum walks. 

An outline of the remainder of the paper is as follows. In Section [5] we define 
the QRW and summarize some known results. Section [3] is concerned with one- 
dimensional QRWs. We develop some theoretical results specific to one dimension, 
that hold for an arbitrary number of chiralities. We work an example to illustrate 
the new phenomena as well as some techniques of computation. Section |4] is con- 
cerned with examples in two dimensions. In particular, we compute the bounding 
curves for some examples previously examined in [BBBP08) . 



2. Background 

2.1. Construction. To specify a lattice quantum walk one needs the di- 
mension d > 1, the number of chiralities k > d + 1, a sequence of k vectors 
v^'^\ . . . , v^''"' e Z'', and a unitary matrix U of rank k. The state space for the 
QRW is 

n-.^L^ {Z'^ X {!,..., k}) . 

A Hilbert space basis for H. is the set of elementary states 6r.j, as r ranges over 
Z'^ and 1 < j < fc; we will also denote 6rj simply by {r,j). Let I ®U denote 
the unitary operator on U, whose value on the elementary state (r,j) is equal to 
^iji^^ Let T denote the operator whose action on the elementary states is 
given by T{r,j) = (r + v^-^^ j). The QRW operator S = S^j^ u f^^u')} is defined by 

(2.1) 5 :=r •(/»[/). 

2.2. Interpretation. The elementary state (r,j) is interpreted as a particle 
known to be in location r and having chirality j. The chirality is a state that can 
take k values; chirality and location are simultaneously observable. Introduction 
of chirality to the model is necessary for the existence of nontrivial translation- 
invariant unitary operators, as was observed by |Mey96| . A single step of the 
QRW consists of two parts: first, leave the location alone but modify the state 
by applying [/; then leave the state alone and make a determinstic move by an 
increment, v^-') corresponding to the new chirality, j. The QRW is translation 
invariant, meaning that if a is any translation operator (r, j) i— > (r + u,j) then 
S o a = (7 o S . The ri-step operator is 5". Using bracket notation, we denote the 
amplitude for finding the particle in chirality j and location x + r after n steps, 
starting in chirality i and location x, by 

(2.2) a{i,j,n,v) ((x, z) |5"|(x + r, j)) ■ 

By translation invariance, this quantity is independent of x. The squared modulus 
|a(i, j, n, r)p is interpreted as the probability of finding the particle in chirality j and 
location x + r after n steps, starting in chirality i and location x, if a measurement 
is made. Unlike the classical random walk, the quantum random walk can be 
measured only at one time without disturbing the process. We may therefore study 
limit laws for the quantities a{i,j,n,r) but not joint distributions of these. 
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2.3. Generating functions. In what follows, we let x denote the vector 
(xi, . . . ,Xd)- Given a lattice QRW, for I < i, j < k we may define a power se- 
ries in d + 1 variables via 

(2.3) ^y(x,y) ^ ^ a(i,j>,r)xV- 

Here and throughout, x"" denotes the monomial power • • ■ x^/ . We let F denote 
the generating matrix {Fij)i<ci_j<k^ which is a fcx fc matrix with entries in the formal 
power series ring in d + 1 variables. The following result from |BP07) is obtained 
via a straightforward use of the transfer matrix method. 

Lemma 2.1 f |BP07[ Proposition 3.1]). Let Af(x) denote the kxk diagonal matrix 
whose diagonal entries are x'^* \ ■ ■ ■ , x"^* ' . Then 



(2.4) F{^,y)^iI-yM{x.)Uy 
Consequently, there are polynomials Pij (x, y) such that 

(2.5) F., = ^ 
where Q{yi, y) := dct(/ - yM{x)U). 

□ 

Let z denote the vector (x, y) and let 

V := {z e C^+i : Q(z) - 0} 

denote the algebraic variety which is the common pole of the generating functions 
Fij. Let Vi := V n T'^^^ denote the intersection of the singular variety V with the 
unit torus T''^^ :— {\xi\ = ■ ■ ■ = \xd\ = \y\ — 1}. An important map on V is the 
logarithmic Gauss map /i : V ^ CP'' defined by 

(2.6) /i(z) := zi— : ... : Zd+i- 



' dzi dzd+i 

The map ^ is defined only at points of V where the gradient VQ does not vanish. 
In this paper we will be concerned only with instances of QRW satisfying 

(2.7) VQ vanishes nowhere on Vi . 

This condition holds generically. 

2.4. Knov^rn results. It is shown in [BBBPOSl Proposition 2.1] that the 
image /i[Vi] is contained in the real subspace MP'' C CP''. Also, under the hypoth- 
esis (|2.7p . dQ/dy cannot vanish on Vi, hence we may interpret the range of /i as 

C RP'' via the identification (xi : ■ ■ ■ : Xd : y) ^ {{xi/y), . . . , (xd/y))- In what 
follows, we draw heavily on two results from [B BBP08| . 

Theorem 2.2 (shape theorem |BBBP08[ Theorem 4.2]). A ssume (j2.7p and let 

G C R'' he the closure of the image of fi on Vi. If K is any compact subset ofG'^, 
then 

a(z,j,n,r) = 0(e--") 
for some c — c{K) > 0, uniformly as r/n varies over K . 
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□ 

In other words, under ballistic rescaling, the region of non-exponential decay 
or feasible region is contained in G. The converse, and much more, is provided by 
the second result, also from the same theorem. For z € Vi, let k(z) denote the 
curvature of the real hypersurface — ilogVi C Mf^^^ at the point logz, where log is 
applied to vectors coordinatewise and manifolds pointwise. 



Theorem 2.3 (asymptotics in the feasible region). Suppose Q satisfies p.7p . For 
r e G, let Z{r) denote the set fj,^^{r) of pre-images in Vi of the projective point r 
under fj,. If k[z) ^ for all z G Z{r), then 



(2.8) a(i,j,n,r) = n-'^/2 



E 

zeZ{r) 



|ViogQ(z)| 



|k(z)| 



-l/2g2cj(r,7i) 



where the argument w(r, n) is given by — r • Arg (z) + i7rr(z)/4 and t{z) is the index 
of the quadratic form defining the curvature at the point (l/i) logz e (l/i) log Vi. 



□ 



3. One-dimensional QRW with three or more chiralities 



3.1. Hadamard QRW. The Hadamard QRW is the one-dimensional QRW 
with two chiralities that is defined in [ADZ93] and analyzed in [ABN+Ol] and [CIR03] . 



It has unitary matrix U = —= 



1 

3 



which is a constant multiple of a 



1 1 

1 -1 

Hadamard matrix, these being matrices whose entries are all ±1. Applying an 
affine map to the state space, we may assume without loss of generality that the 
steps are and 1. Up to a rapidly oscillating factor due to a phase difference in 
two summands in the amplitude, it is shown in these early works that the rescaled 
amplitudes n^^'^aii, j,n,n6) converge to a profile f{9) supported on the interval 
1 \/2 1 V2 



J := 



4 ' 2 



The function / is continuous on the interior of J and 



blows up like \9 — 6q\~'^^'^ when 9q is an endpoint of J. These results are extended 
in |BP07j to arbitrary unitary matrices. The limiting profiles are all qualitatively 
similar; a plot for the Hadamard QRW is shown in figure[Tl with the upper envelope 
showing what happens when the phases of the summands line up. 



3.2. Experimental data with three or more chiralities. When the num- 
ber of chiralities is allowed to exceed two, new phenomena emerge. The possibility 
of a bound state arises. This means that for some fixed location a;, the ampli- 
tude a{i,j,n,x) does not go to zero as n ^ oo. This was first shown to occur 
in [BCA03"l IKonOS] . From a generating function viewpoint, bound states occur 
when the denominator Q of the generating function factors. The occurrence of 
bound states appears to be a non-generic phenomenon. 

In 2007, two freshman undergraduates, Torin Greenwood and Rajarshi Das, 
investigated one-dimensional quantum walks with three and four chiralities and 
more general choices of U and {v*-^-*}. Their empirical findings are catalogued at 

\protect\vrule widtliOpt\protect\href {http : //www . math . upenn . edu/\string~pemantle/Summer2007/Fi 
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Figure 1. probability profile for the 1-dimensional Hadamard 
QRW: red line is upper envelope, computed by aligning the phases 
of the summands, while black points are actual squared magnitudes 



The probability profile shown in figure [5] is typical of what they found and is the 
basis for an example running throughout this section. In this example, 

17 6 20 -2 ' 
-20 12 13 -12 
-2 -15 4 -22 
-6 -18 12 15 

and v^-') = —1, 0, 1, 2 for j — 1, 2, 3, 4 respectively. The profile shown in the figure is 
a plot of |a(l, 1, 1000, x)\'^ against x for integers x in the interval [—1000, 2000]. The 
values were computed exactly by recursion and then plotted. The most obvious new 
feature is the existence of a number of peaks in the interior of the feasible region. 
The phase factor is somewhat more chaotic as well, which turns out to be due to 
a greater number of summands in the amplitude function. Our aim is to use the 
theory described in Section [2] to establish the locations of these peaks, that is to 
say, the values of 6 for which n^/'^a{i, j,n, x) become unbounded for x sufficiently 
near n9. 

3.3. Results and conjectures. The results of Section[2]may be summarized 
informally in the case of one-dimensional QRW as follows. Provided the quantities 
VQ and k do not vanish for the points z associated with a direction r, then the 
amplitude profile will be a the sum of terms whose phase factors may be somewhat 
chaotic, but whose magnitudes are proportional to ViogQ|. In practice the 

magnitude of the amplitude will vary between zero and the sum of the magnitudes 
of the pieces, depending on the behavior of the phase terms. In the two-chirality 
case, with only two summands, it is easy to identify the picture with the theoretical 



(3.1) 



U = 



27 
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Figure 2 . probability profile for a four-chirality QRW in one dimension 

result. However, in the multi-chirality case, the empirical results in figureOare not 
easily rectified with the theoretical result, firstly because the theoretical result is 
not trivial to compute, and secondly because the computation appears at first to be 
at odds with the data. In the remainder of Section [Sj we show how the theoretical 
computations may be executed in a computer algebra system, and then recify these 
with the data in figure [21 The first step is to verify some of the hypotheses of 
Theorems 12.2112.31 The second step, reconciling the theory and the data, will be 
done in Section [3^ 

Proposition 3.1. Let Q(x, y) he the denominator of the generating function for 
any QRW in any dimension that satisfies the smoothness hypothesis (|2.7p . Let tt 
be the projection from Vi to the d-torus T'^ that forgets the last coordinate. Then 
the following properties hold. 

(i) dQ/dy does not vanish on Vi; 

(ii) Vi is a compact d-manifold; 

(iii) TT : Vi — !■ T'^ is smooth and nonsingular; 

(iv) In fact, Vi is homeomorphic to a union of some number s of d-tori, each 
mapping smoothly to T'^ under tt and covering T'^ some number nj times 
for 1 < J < s. 

(v) K : Vi ^ M vanishes exactly when the determinant of the Jacobian of the 
map fj, vanishes. 

(vi) K vanishes on the boundary dfi[Vi] of the range of fi. 

Proof. The first two conclusions are shown as |BBBP08| Proposition 2.2]. 
The map tt is smooth on T'^'^^, hence on Vi, and nonsingularity follows from the 
nonvanishing of the partial derivative with respect to y. The fourth conclusion 
follows from the classification of compact d-manifolds covering the d-torus. For the 
fifth conclusion, recall that the Gauss-Kronecker curvature of a real hypersurface is 
defined as the determinant of the Jacobian of the map taking p to the unit normal 
at p. We have identified projective space with the slice Zd+i = 1 rather than with 
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the slice |z| ~ 1, but these are locally diffeomorphic, so the Jacobian of /i still 
vanishes exactly when k vanishes. Finally, if an interior point of a manifold maps 
to a boundary point of the image of the manifold under a smooth map, then the 
Jacobian vanishes there, hence the last conclusion follows from the fifth. □ 

An empirical fact is that in all of the several dozen quantum random walks we 
have investigated, the number of components of Vi and the degrees of the map tt 
on each component depend on the dimension d and the vector of chiralitics, but 
not on the unitary matrix U . 

Conjecture 3.2. If d, fc, v'^^\ . . . , v^'^^ are fixed and U varies over unitary matri- 
ces, then the number of components of V and the degrees of the map tt on each 
component are constant, except for a set of matrices of positive co-dimension. 

Remark. The unitary group is connected, so if the conjecture fails then a transition 
occurs at which Vi is not smooth. We know that this happens, resulting in a bound 
state |Kon05| . however in the three-chirality case, the degeneracy does not seem 
to mark a transition in the topology of Vi . 

Specializing to one dimension, the manifold Vi is a union of topological circles. 
The map /i : Vi — > M is evidently smooth, so it maps Vi to a union of intervals. In 
all catalogued cases, in fact the range of ^ is an interval, so we have the following 
open question: 

Question 3.3. Is it possible for the image of /i to be disconnected? 

Because /i smoothly maps a union of circles to the real line, the Jacobian of the 
map /i must vanish at least twice on each circle. Let W denote the set of z G Vi for 
which k(z) = 0. The cardinality of W is not an invariant (compare, for example, 
the example in Section 13.41 with the first 4-chirality example on the web archive). 
This has the following interesting consequence. Again, because the unitary group 
Uk is connected, by interpolation there must be some U for which there is a double 
degeneracy in the Jacobian of /i. This means that the Taylor series for logy on Vi 
as a function of log x is missing not only its quadratic term but its cubic term as 
well. In a scaling window of size v}!"^ near the peaks, it is shown in [BP07| that the 
amplitudes are asymptotic to an Airy function. However, with a double degeneracy, 
the same method shows a quartic-Airy limit instead of the usual cubic- Airy limit. 
This may be the first combinatorial example of such a limit and will be discussed 
in forthcoming work. 

Let W = {w*^°\ . . . ,w(*)} be a set of vectors in R". Say that W is rationally 
degenerate if the set of i-tuples (r • (w — w(''))wgM/ is not dense in (R mod 27r)* as 
r varies over Z" . Generic i-tuples are rationally nondegenerate because degeneracy 
requires a number of linear relations to hold over the 27r(Q. If W is rationally nonde- 
generate, then the distribution on t-tuples (r- (w — w('^^)weH' when r is distributed 
uniformly over any cube of side M in Z'* converges weakly to the uniform distri- 
bution on (Z mod 27r)*. Let x(ai,...,at) denote the distribution of the square 
modulus of the sum of t complex numbers chosen independently at random with 
moduli ai,...,Q!t and arguments uniform on [— 7r,7r]. The following result now 
follows from the above discussion. Theorems 12.21 and 12.31 and Proposition 13. II 

Proposition 3.4. For any one- dimensional QRW, let Q, Z{r) and k be as above. 
Let J be the image ofVi under /i. Let r be any point of J such that k{z) ^ for all 
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z G Z{r) and W :— (1/i) logZ(r) is rationally nondegenerate. Then for any e > 
there exists an M such that ifrin) is a sequence of integer vectors with r*^") /n — > r, 
the empirical distribution of n'^ times the squared moduli of the amplitudes 

{a{i,j, n, r{n) + e) : C e {0, . . . , M - 1}^+'} 

is within e of the distribution x(q^i, • ■ ■ , Q^t) where t = |Z(r)|, {z*^^'} enumerates 
Z(r), and aj — |_Py (z'^-')k(z^^-')~"'^/^|. // r ^ J, then the empirical distribution 
converges to a point mass at zero. 

□ 

Remark. Rational nondegeneracy becomes more difRcult to check when the size 
of Z{y) increases, which happens when the number of chiraUties increases. If one 
weakens the conclusion to convergence to some nondegenerate distribution with 
support in / := [0, ^ |Py(z)^K(z)~-'^|], then one needs only that not all components 
of all differences logz — logz' are rational, for z,z' e Z{y). For the purpose of 
qualitatively explaining the plots, this is good enough, though the upper envelope 
may be strictly less than the upper endpoint of / (and the lower envelope may be 
strictly greater than zero) if there is rational degeneracy. 

Comparing to figure [21 we see that J appears to be a proper subinterval of 
[—1,2], that there appears to be up to seven peaks which are local maxima of 
the probability profile. These include the endpoints of J (cf. the last conclusion 
of Proposition 13. ip as well as several interior points, which we now understand 
to be places where the map pL folds back on itself. We now turn our attention 
to corroborating our understanding of the picture by computing the number and 
locations of the peaks. 

3.4. Computations. Much of our computation is carried out symbolically in 
Maple. Symbolic computation is significantly faster when the entries of U are ratio- 
nal, than when they are, say, quadratic algebraic numbers. Also, Maple sometimes 
incorrectly simplifies or fails to simplify expressions involving radicals. It is easy 
to generate quadratically algebraic orthogonal or unitary matrices via the Gram- 
Schmidt procedure. For rational matrices, however, we turn to a result we found 
in [L091j . 

Proposition 3.5. The map S (I + S){I — S)^^ takes the skew symmetric 
matrices over a field to the orthogonal matrices over the same field. To generate 
unitary matrices instead, use skew-hermitian matrices S. 

□ 

The map in the proposition is rational, so choosing S to be rational, we obtain 
orthogonal matrices with rational entries. In our running example, 

"0 -3-1 3 " 

3 1-2 

1-10 2 

-3 2-2 

leading to the matrix U of equation (|3.ip . 
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The example shows amplitudes for the transition from chirality 1 to chirality 1, 
so we need the polynomials Pn and Q: 

Pii{x,y) = (27a; - 15ya;^ - 4?;a; + 12?/V - 12 2; + 42;V + 9?;2 - ITyV) a: 
Q{x,y) = -ITyV + 9t/2 + 27a;- 12y + 12y2a;3 + 8y2a;2 - ISyx^ - 4y^a;^ 
-15 y^x + 12 y^x - Ayx - 17 yx^ + 9 y^x* - 12 y^x^ + 27 y^x^ . 
The curvature is proportional to the quantity 

( a; Qx y Qy) Qx y Qy y {Qy Qxx ~t~ Qx Qxy 2 Qx Qy Qxy) ; 

where subscripts denote partial derivatives. Evaluating this leads to xy times a 
polynomial K{x,y) that is about half a page in Maple 11. The command 

Basis ( [Q, if] , plex (y,x)); 
leads to a Grobner basis, the first element of which is an elimination polynomial 
p{x), vanishing at precisely those x-values for which there is a pair {x,y) G V for 
which K{x,y) — 0. We may also verify that Q is smooth by computing that the 
ideal generated by [Q, Qx, Qy] has the trivial basis, [1]. 

To pass to the subset of roots of p{x) that are on the unit circle, one trick is 
as follows. If z = X + 1/x then x is on the unit circle if and only if z is in the 
real interval [—2,2]. The polynomial defining z is the elimination polyomial q{z) 
for the basis [p, 1 ~ zx + x^\. Applying Maple's built-in Sturm sequence evaluator 
to q shows symbolically that there are six roots of z in [—2,2]. This leads to six 
conjugate pairs of x values. The second Grobner basis element is a polynomial 
linear in j/, so each x value has precisely one corresponding y value. The y value 
for X is the conjugate of the y value for a;, and the function [i, takes the same value 
at both points of a conjugate pair. Evaluating the [i function at all six places leads 
to floating point expressions approximately equal to 

1.362766, 1.126013, 0.929248, 0.229537, -0.143835, -0.346306 . 

Drawing vertical lines corresponding to these six peak locations leads to figure [3l 
Surprisingly, the largest peak appearing in the data plot appears to be missing 
from the set of analytically computed peak directions. Simultaneously, some of the 
analytically computed peaks appear quite small and it seems implausible that the 
probability profile blows up there. Indeed, this had us puzzled for quite a while. 
In order to doublecheck our work, we plotted y against a;, resulting in the plot in 



figure 4(a) which should be interpreted as having periodic boundary conditions 
because x and y range over a circle. This shows Vi to be the union of two circles, 
each embedded in so that the projection tt onto x has degree 2. (Note: the 
projection onto y has degree 1, and the homology class of the embedded circle is 
(2,-1) in the basis generated by the x and y axes.) We also plotted [i against 
X. To facilitate computation, we used Grobner bases to eliminate y from Q and 
xQx — ^^yQy^ enabling us to plot solutions to a single polynomial. The resulting 
plot is shown in figure [4(b)| 

The last figure shows nicely how peaks occur at values where the map /i back- 
tracks. The explanation of the appearance of the extra peak at « 0.7 becomes 
clear if we compare plots at n = 1,000 and n — 10,000. At first glance, it looks 
as if the extra peak is still quite prominent, but in fact it has lowered with respect 
to the others. To be precise, the false peak has gone down by a factor of 10, from 
0.004 to 0.0004, because its probabilities scaled as n~^. The width of the peak also 
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Time 



Figure 3. probability profile with peaks drawn as vertical lines 





(a) y versus x 



(b) ^ versus x 



Figure 4. Two interleaved circles and their images under the 
Gauss map 



remained the same, indicating convergence to a finite probability profile. The real 
peaks, however, have gone down by factors of lO^^'^, as is shown to occur in the 
Airy scaling windows near directions r where k{z) = for some z G Z{y). When 
the plot is vertically scaled so that the highest peak occurs at the same height in 
each picture, the width above half the maximum has shrunk somewhat, as must 
occur in an Airy scaling window, which has width ^Jn. The location of the false 
peak is marked by a nearly flat spot in figure [4(b)| at height around 0.7. The curve 
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2,000 -10,000 




(a) n = 1000 (b) n = 10000 

Figure 5. As n — > cxd, one peak scales down more rapidly 



stays nearly horizontal for some time, causing the false peak to remain spread over 
a macroscopic rescaled region. 



4. Two-dimensional QRW 

In this section we consider two examples of QRW with d = 2, fc = 4 and 
steps v(i) = (0,0), v(2) = (1,0), = (0,1) and = (l,l). To complete the 
specification of the two examples, we give the two unitary matrices: 

1 1 11' 

-11-11 
1 -1-11 
-1-111 

1 1 1 1 " 

-11-11 
-11 1 -1 ■ 
-1-11 1 

Note that these are both Hadamard matrices; neither is the Hadamard matrix with 
the bound state considered in [M00O4] , nor is either in the two-parameter family 
referred to as Grover walks in [WKKKOSj . The second differs from the first in that 
the signs in the third row are reversed. Both are members of one-parameter families 
analyzed in |BBBP08] . in Sections 4.1 and 4.3 respectively. The (arbitrary) names 
given to these matrices in jBra07|, [BBBP08 are respectively 5(1/2) and i3(l/2). 
Intensity plots at time 200 for these two quantum walks, given in figure[Sl reproduce 
those taken from [BBBPOS] but with different parameter values (1/2 each time, 
instead of 1/8 and 2/3 respectively). 

For the case of C/i it is shown in [BBBPOSl Lemma 4.3] that Vi is smooth. 
Asymptotics follow, as in Theorem 12.31 of the present paper, and an intensity plot 
of the asymptotics is generated that matches the empirical time 200 plot quite well. 



(4.1) 



- 2 



(4.2) 



U2 
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(a) Ui 



(b) U2 



Figure 6. Time 200 probability profiles for two quantum walks: 
the darkness at (r, s) corresponds to the squared amplitude 
|a(l,l,200,r,s)|2. 



In the case of U2, Vi is not smooth but jBBBPOSl Theorem 3.5] shows that the 
singular points do not contribute to the asymptotics. Again, a limiting intensity 
plot follows from Theorem l2.3l of the present paper and matches the time 200 profile 
quite well. 

It follows from Proposition 13.41 that the union of darkened curves where the 
intensity blows up is the algebraic curve where k vanishes, and that this includes the 
boundary of the feasible region. The main result of this section is the identification 
of the algebraic curve. While this result is only computational, it is one of the 
first examples of computation of such a curve, the only similar prior example being 
the computation of the "Octic circle" boundary of the feasible region for so-called 
diabolo tilings, identified without proof by Cohn and Pemantle and first proved 
by [KO07| (see also |BP10| ). The perhaps somewhat comical statement of the 
result is as follows. 

Theorem 4.1. For the quantum walk with unitary coin flip U2, the curvature of 
the variety Vi vanishes at some z e Z(r, s) if and only if (r, s) is a zero of the 
polynomial P2 and satisfies \r\ + \s\ < 3/4, where 

P^(r, s) := 1 + 14r2 - 3126r4 + 97752r6 - 1445289r* + 12200622ri" - 6Al50356r^^ + 
220161216ri4 - 504431361ri6 + 774608490ri8 - 785130582r20 + 502978728r22 - 
184298359r24+29412250r26+14s2-1284r2s2_ii30l6r4s2+5220612rS2-96417162r8s2 + 

924427224riOs2-4865103360ri2s2+14947388808ri''s2_27714317286ri6s2+30923414124ri8s2- 
19802256648r20s2 4. 6399721524r22s2 _ 721963550r24s2 _ 31265* - 113016r2s4 + 
7942218r4s^-68684580rS^-666538860rSs*+15034322304ri°s''-86727881244ri2s4+ 
226469888328ri4s4 - 296573996958riS4 + 183616180440ri*^s4 _ 32546593518r2"s4 _ 
8997506820r22s4+97752s64_5220612r2s6_68684580r4s6+3243820496rSS-25244548160r«sS + 
59768577720riOs6 - 147067477144ri2s6 +458758743568ri4s6 - 749675452344riS6 + 
435217945700ri8s6-16479111716r2"s'5-1445289s*^-96417162r2s8-666538860r4s*- 
25244548160r*^s^ + 194515866042r«s8 - 421026680628ri0s8 + 611623295476ri2s8 - 
331561483632ri4s8+7820601831riS*^+72391117294ri8s*+12200622si°+924427224r2.si"+ 
15034322304r4sio + 59768577720rSi'' - 421026680628r**si° + 421043188488riOsi° - 
1131276050256ri2sio - 196657371288ri4sio + 151002519894riSi° - 64150356si2 _ 
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4865103360r2si2 - 86727881244r4si2 _ l47067477144rSi2 + 611623295476r8si2 _ 
1131276050256ri''si2 + 586397171964ri2si2 _ 231584205720ri4si2 -|- 220161216si'' + 
14947388808r2s" + 226469888328r4si4 +458758743568r6si4 - 331561483632r«si4 - 
196657371288riOsi'* - 231584205720ri2si4 _ 504431361si6 - 27714317286r2si6 - 
296573996958r4si6 _ 749675452344r6si6 + 7820601831r*^si6 + 151002519894ri°si6 + 
774608490si8+30923414124r2si8_|_i836i6i80440rSi8+435217945700rSi8+72391117294rS 
785130582s20_i9802256648r2s20_32546593518r4s20-16479111716r6s20+502978728s22+ 
6399721524r2s22-8997506820r4s22-184298359s24-721963550r2s24+29412250s26. 

We check visually that the zero set of P2 does indeed coincide with the curves 
of peak intensity for the U2 QRW. 




(a) probabilities at time 200 at (r, s) (b) zero set of P2 in {r/n,s/n) 

Figure 7. The probability profile for the U2 QRW alongside the 
graph of the zero set of P2 

Before embarking on the proof, let us be clear about what is requred. If r 
is in the boundary of the feasible region, then k must vanish at the pre-images 
of r in the unit torus. The boundary, dG, of the feasible region is therefore a 
component of a real algebraic variety, W. The variety W is the image under the 
logarithmic Gauss map fi of the points of the unit torus where Q and k both 
vanish. Computing this variety is easy in principle: two algebraic equations in 
(x, y, z, r, s) give the conditions for fj,{x, y, z) — (r, s) and two more give conditions 
for Q{x, y, z) = k(x, y, z) — 0; algebraically eliminating {a:, y, z} then gives the 
defining polynomial P2 for W . In fact, due to the number of variables and the 
degree of the polynomials, a straightforward Grobner basis computation does not 
work and we need to use iterated resultants in order to get the computation to halt. 
The last step is to discard extraneous real zeros of P2 , namely those in the interior 
of G or G^, so as to arrive at a precise description of dG. 
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Proof. To eliminate subscripts, we use the variables (x, y, z) instead of (xi, X2, y)- 
The condition for z G Z{r, s) is given by the vanishing of two polynomials Hi and 
H2 in {x,y,z,r,s), where 

Hi{x,y,z,r,s) := xQ^-rzQ^; 
H2{x,y,z,r,s) := yQy-szQ^. 

The curvature of Vi at z also vanishes when a single polynomial vanishes, which 
we will call L(x,y,z). While explicit formulae for L may be well known in some 
circles, we include a brief derivation. For {x,y,z) £ Vi, write x = e*^,y = 
e'^ and z — e*'^ . By Proposition 13.11 we know that Qz ^ on Vi , hence the 
parametrization of Vi by X and Y near a point {x,y,z) is smooth and the par- 
tial derivatives Zx, Zy, Zxx, Zxy Zyy are well defined. Implicitly differentiating 
(5(e*^, e''^, e'^t^'"*')) = with respect to X and Y we obtain 

^^ = -^ and = 
and differentiating again yields 



-ixz 



^xx 



jyy — 



{zQzf 

-iyz 



[QxQzizQz - 'ixzQxz + xQx) + Xz{QIQzz + QIQxx)] 
3 [QyQz{zQz - '^yzQyz + zQy) + yz{QyQzz + QlQyy)] 



{zQz 

'ixyz 

In any dimension, the Gaussian curvature vanishes exactly when the determinant 
of the Hessian vanishes of any parametrization of the surface as a graph over d — 1 
variables. In particular, the curvature vanishes when 

Zxx Zxy 
Zxy Zyy 

vanishes, and plugging in the computed values yields the polynomial 



det 



L{x,y,z) := -xyzQlQly + zQxQlQy-2yzQxQzQyQyz+yQxQzQl + yzQxQlQzz 

-\- y zQ xQ zQyy ^"^zQxQ zQxzQy ~l~ ^-^yzQxQxzQyQyz '^^yzQxQzQxzQyy 
~\~ xQ^QzQy H" "^yQxQzQyy -^^Q xQ zzQy H~ -^V^QxQzzQyy XzQxxQ zQy 

- 2xyzQ 

xxQ zQyQyz H" ^yQxxQzQy H" -^yzQxxQyQ zz xy zQ xxQ zQ yy 

-xyzQyQ^^ 

- XyzQ^Qy^ + 2xyzQ zQxyQxQyz + '^XyzQ zQxyQyQxz ' - 2xyQzQ xyQxQy 

- 2xyzQ 

xyQxQyQ zz ' 

It follows that the curvature of Vi vanishes for some {x,y,z) G Z{r,s) if and 
only if the four polynomials Q, Hi, H2 and L all vanish at some point {x, y, z, r, s) 
with {x,y, z) G . I gnoring the condition {x,y,z) G for the moment, we see 
that we need to eliminate the variables (x, y, z) from the four equations, leading 
to a one-dimensional ideal in r and s. Unfortunately Grobner basis computations 
can have very long run times, with published examples showing for example that 
the number of steps can be doubly exponential in the number of variables. Indeed, 
we were unable to get Maple to halt on this computation (indeed, on much smaller 
computations). The method of resultants, however, led to a quicker elimination 
computation. 
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Definition 4.2 (resultant). Let f{x) :— Y^j^o^i-'^'' ^^"^ di^) '■~ Sj=o ^i^"' 
two polyomials in the single variable x, with coefficients in a field K. Define the 
resultant firesult(/, g, x) to be the determinant of the {£ + rn) x {£ + m) matrix 



/ ao 

ai 


ao 




bo 
bi 


bo 




a2 


ai 




62 


bi 






a2 


■ ao 




b2 ' 


■ bo 


ai 




■ ai 






■ 61 




ai 


a2 






b2 


\ 




ai 






bra 



The crucial fact about resultants is the following fact, whose proof may be 
found in a number of places such as |CL098|, [GKZ94) : 

(4.3) 6result(/, g,x)^0 ^ 3x : f{x) ^ g{x) ^ . 

Iterated resultants are not quite as nice. For example, if /, h are polynomials in 
X and ?/, they may be viewed as polynomials in y with coefficients in the field of 
rational functions, K{x). Then firesult(/, /i, y) and firesult(g, /i, y) are polynomials 
in X, vanishing respectively when the pairs (/, h) and (5, h) have common roots. 
The quantity 

R fSresult(6result(/, /i, y), f5result((;, h, y), x) 

will then vanish if and only if there is a value of x for which /(a;, yi) = h{x, yi) — 
and g{x, 1/2) — h{x, ^2) = 0. It follows that if f{x, y) — g{x, y) = then R = 0, but 
the converse does not in general hold. A detailed discussion of this may be found 
in [BM07I . 

For our purposes, it will sufhce to compute iterated resultants and then pass 
to a subvariety where a common root indeed occurs. We may eliminate repeated 
factors as we go along. Accordingly, we compute 

R12 ;= TZad{&resu\t{Q,L,x)) 
Ri3 ; = 7^ac^(firesuh (0 , Hi,x)) 
Rii ;= 7^ac^(firesuh(Q,i^2,a;)) 

where TZad{P) denotes the product of the first powers of each irreducible factor of 
P. Maple is kind to us because we have used the shortest of the four polynomials, 
Q, in each of the three first-level resultants. Next, we eliminate y via 

i?i24 '■— 7^ad(firesult(i?i2, -R14, y)) 
i?i34 :— 7^ad(firesult(i?i3, i?i4, y)) . 

Polynomials i?i24 and i?i34 each have several small univariate factors, as well as 
one large multivariate factor which is irreducible over the rationals. Denote the 
large factors by /124 and /134. Clearly the univariate factors do not contribute to 
the set we are looking for, so we eliminate z be defining 

-R1234 := 7^a(i(6result(/i24, /134, z)) . 
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Maple halts, and we obtain a single polynomial in the variables (r, s) whose zero set 
contains the set we are after. Let Q denote the set of (r, s) such that k{x, y, z) = 
for some {x,y,z) e V with fi{x,y,z) = (r, s) [note: this definition uses V instead 
of Vi.] It follows from the symmetries of the problem that il is symmetric under 
r I— > — r as well as s i— s- — s and the interchange of r and s. Computing iterated 
resultants, as we have observed, leads to a large zero set ft'; the set 57' may not 
possess r-s symmetry, as this is broken by the choice of order of iteration. Factoring 
the iterated resultant, we may eliminate any component of Q' whose image under 
transposition of r and s is not in Q' . Doing so, yields the irreducible polynomial 
P2- Because the set is algebraic and known to be a subset of the zero set of the 
irreducible polynomial P2, we see that ft is equal to the zero set of P2. 

Let rip ^ denote the subset of those (r, s) for which as least one {x, y, z) G 
fj,^^{{r, s)) with K(x,y,z) = lies on the unit torus. It remains to check that flo 
consists of those (r, s) e with \r\ + \s\ < 3/4. 

The locus of points in V at which k vanishes is a complex algebraic curve 7 
given by the simultaneous vanishing of Q and L. It is nonsingular as long as VQ 
and VL are not parallel, in which case its tangent vector is parallel to VQ x VL. 
Let p := xQx/{zQz) and a := yQy/{zQz) be the coordinates of the map /i under 
the identification of CP^ with {(r, s, 1) : r, s e C}. The image of 7 under /i (and this 
identification) is a nonsingular curve in the plane, provided that 7 is nonsingular 
and either dp or da is nonvanishing on the tangent. For this it is sufficient that one 
of the two determinants det Mp, det does not vanish, where the columns of Mp 
are VQ, VL, Vp and the columns of are VQ, VL, Vcr. 

Let (xq, yo, zq) be any point in Vi at which one of these two determinants does 
not vanish. It is shown in [BBBPOS, Proposition 2.1] that the tangent vector to 
7 at {xo,yo,zo) in logarithmic coordinates is real; therefore the image of 7 near 
(xQ,yo, Zq) is a nonsingular real curve. Removing singular points from the zero 
set of P2 leaves a union U of connected components, each of which therefore lies 
in ilo or is disjoint from CIq. The proof of the theorem is now reduced to listing 
the components, checking that none crosses the boundary \r\ + \s\ — 3/4, and 
checking Z(r, s) for a single point (r, s) on each component (note: any component 
intersecting {\r\ + \s\ > 1} need not be checked as we know the coefficients to be 
identically zero here). □ 

We close by stating a result for Ui , analogous to Theorem 14.11 The proof is 
entirely analogous as well and will be omitted. 

Theorem 4.3. For the quantum walk with unitary coin flip Ui, the curvature of 
the variety Vi vanishes at some (x, y, z) € Z{r^ s) if and only if \r\ and \s\ are both 
at most 2/3 and (r, s) is a zero of the polynomial 

Pi(r,s) := 132019ri6 + 2763072s VO - 513216s - 6505200sVi8 + 2bQs\'^ + 
8790436sVi6 10639416siV + 39759700siV'' - 12711677siV + AlA02b7s^\'^ - 
513216s22r2-7492584s2ri4+2503464si'^r6-62208s22+16s6+141048r20+8790436si6r2 + 
2763072s2"r2 - 6505200siV - 40374720siV + 64689624siV - 33614784siV + 
14725472siOri0+121508208si6r«-1543si''-23060sV + 100227200si°ri2+7363872s2M- 
176524ri8 + 121508208s**ri6 - m2nbb2s^r^^ - 13374107sV + 1647627sV + 
18664050s V - 227481984siOri4 - 19343sV + 2792MA%Qs^\^'^ - 67173440si V - 
7492584siV+4140257s2ri2+291173sV-1449662sVio+7363872sV"-227481984siVi" + 
132019si6-197271552siV-59209ri4-1449662siV + 100227200siViO-1543ri"- 
153035200siV-13374107sV+3183044sV+39759700sVi2-176524sis+72718sS-4+ 
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1647627sV - 62208r22 + 141048.s2" - 1472sV + 11664s24 - 33614784.sVi8 + 
128187648siV-1472.s2r4-67173440s4ri4+291173.sV+64689624.s4ri«-10639416s^ 
5920931"* + 72718s''r'5 + 92321584s^ri2 - 56r^ + 92321584si2r« - 153035200s6r" - 
23060sV+128187648sVi'^-40374720s6ri8+72282208siV+14793ri2+11664r24+ 
14793s^2 ^ ^ 2503464s6rio - 56s^ - 12711677s^ri° + 72282208sV^2 

□ 

Summary. We have stated an asymptotic amplitude theorem for general one- 
dimensional quantum walk with an arbitrary number of chiralities and shown how 
the theoretical result corresponds, not always in an obvious way, to data generated 
at times of order several hundred to several thousand. We have stated a general 
shape theorem for two-dimensional quantum walks. The boundary is a part of 
an algebraic curve, and we have shown how this curve may be computed, both in 
principle and in a Maple computation that halts before running out of memory. 
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